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Abstract. Ulysses measurements have confirmed that there is 
no significant gradient with respect to hcliomagnetic latitude in 
the radial component, B rt of the interplanetary magnetic field. 
There are two processes responsible for this observation. In the 
corona, the plasma fi is 1, except directly above streamers, 
so both longitudinal and latitudinal (meridional) gradients in 
field strength will relax, due to the transverse magnetic pressure 
gradient force, as the solar wind carries magnetic flux away from 
the Sun. This happens so quickly that the field is essentially 
uniform by 5 R ( :y. Beyond 10 /?(.,, fi > 1 and it is possible for 
a meridional thermal pressure gradient to redistribute magnetic 
flux - an effect apparently absent in Ulysses and earlier ICE and 
IMP data. Wc discuss this second effect here, showing that its 
absence is mainly due to the perpendicular part ol the anisotropic 
thermal pressure gradient in the interplanetary medium being 
loo small to drive significant meridional transport between the 
Sun and ~ 4 AU. This is done using a linear analytic estimate of 
meridional transport. The first effect was discussed in an earlier 
paper. 
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1. Introduction 

Ulysses in 1993-1995 (Smith and Balogh, 1995; Balogh el al. 
1995) and ICE and IMP-8 in 1984-1988 (Burton et al. 1995) 
observed no significant gradient in hcliomagnetic latitude in the 
radial component of the interplanetary magnetic field (IMF) at 
heliocentric distances of 1 -4 AU. These observations were made 
near solar minimum when the Sun’s magnetic field is nearly an 
axially aligned dipole. Fig. 1 shows the results from IMP-8 
and ICE scaled to 1 AU and plotted against magnetic latitude. 
The data arc five degree latitude bin averages. The solid line 
is a simple equatorial current sheet model of the IMF which, 
outside a few tens of solar radii, gives a uniform field in the two 


hemispheres (Wolfson, 1985) and is a good approximation to 
the radial component, B ry of the IMF at those distances. The 
apparent appearance of a small gradient in the radial component 
of the field near the equator in Fig. 1 is probably all due to small 
errors in magnetic sector identification so that the true gradient 
is completely negligible. No gradient at all was seen by Ulysses 
in 1993-1995 up to heliographic latitudes of ±80.22°. 

It was shown in an earlier paper (Suess and Smith, 1996) 
that the field strength at the lop of the corona, ~ 10 7?.q, has 
no significant gradient in latitude or longitude outside of the 
heliospheric plasma sheet (HPS) which exists above streamers 
(Gosling et al. 1981). 

This is a consequence of the plasma fi, the ratio of thermal 
pressure to magnetic pressure, being <C 1 inside 10 B.& except in 
the HPS. Beyond this distance,/^ > 1 , at least to well beyond the 
orbit of Jupiter, except in some transients. Therefore, between 
the top of the corona and several AU, meridional gradients in 
the thermal pressure may redistribute magnetic llux in latitude. 
A linear analytic model is used here to analyze the Ulysses data 
for the effect of meridional gradients in the density, temperature, 
and flow speed on the meridional distribution of magnetic flux. 

It has already been suggested by Suess et al. (1977) that the 
thermal pressure gradient may indeed have a measurable effect. 
Suess et al. utilized a semi-analytic MHD model of coronal hole 
flow to simulate the northern polar coronal hole observed from 
Skylab in 1973 and described by Munro and Jackson (1977). 
Their model, when carried to I AU, implied that meridional re- 
distribution of mass and magnetic flux between the high corona 
and 1 AU might be a significant process. Results shown in Fig. 
2 from the Suess et al. model illustrate the variation of the radial 
magnetic field with colalitude at 2.0 R&, 5.0 and 1 AU. 
Panel (a) shows what is actually the boundary condition at 2 
a smooth, dipole field. To understand panels (b) and (c), a 
little more must be said about the model. Briefly, it is polytropic 
with an index 7 = 1.05 and is therefore nearly isothermal in 
radius, although the temperature varies strongly with latitude. 
The plasma 3 is less than unity in the corona, but is larger than 
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HELIOMAGNETIC LATITUDE 

Fig. 1. Radial magnetic field values in five degree latitude bins, scaled 
to 1 AU, from ICE and IMP-8 during 1984-1988, bracketing solar 
minimum, and Ulysses during 1993-1994, also at the time of solar 
minimum (Burton el al. 1995). The solid line which is discontinuous 
at the equator is from a simple current sheet model (Wolfson, 1985) at 
a heliocentric distance of more than ca. 50 Rq. 

unity beyond 5-10 R s . The effects both of the relaxation, with 
increasing radius, of meridional gradients in magnetic pressure 
in the corona (Suess and Smith, 1996) and of the interplanetary 
thermal pressure gradient can be seen in this figure. Panel (b) 
shows that the field is almost independent of colatitude at 5 /?© 
due to the relaxation of the meridional magnetic pressure gradi- 
ent in the low 3 corona. Panel (c) shows that a large variation in 
B r is reintroduced between 5 Rq and 1 AU by the interplanetary 
meridional thermal pressure gradient. 

While Suess et al. successfully simulated the observed coro- 
nal hole properties, the IMP and ICE results shown in Fig. 1 
demonstrate that the conditions at 1 AU predicted by their model 
did not exist at the time those observations were made, or at the 
time the more recent Ulysses observations were made (Smith 
and Balogh, 1995; Balogh, et al. 1995). The failure of the model 
to reproduce the Ulysses, ICE, and IMP observations is due both 
to the model pressure and the meridional pressure gradient in 
the interplanetary medium being too large. The large pressure 
is a consequence of using 7 = 1 .05 so that the temperature at 1 
AU is between 0.5 x 1 0 6 and 1 .0 x 1 0 6 K - larger than observed. 
The large pressure gradient is due to the large temperature vari- 
ation across the model coronal hole. When combined with the 
predicted density, this leads to a large meridional pressure gra- 
dient that redistributes mass and momentum flux. In principle, 
that model could be utilized again to analyze Ulysses data and 
understand the absence of interplanetary meridional transport. 
However, the problem is physically and mathematically some- 



Polar Angie 

Fig. 2a-c. The radial magnetic field variation between 0 = 0" (north 
pole) and 0 = 60", at a 2 /?©, b 5 Rq, c 1 AU, from the MHD model 
of the 1973 polar coronal hole of Suess et al. (1977). 

what simpler than modeling coronal flow where radial gradients 
in the flow speed are large. Therefore, we instead describe a sim- 
pler, fully analytic analysis of the data. 

The result we find is suprising. Using the simplest assump- 
tions that the pressure is isotropic and that the meridional pres- 
sure gradient acts uniformly all the way from 10 /?© to 4 AU 
leads to a contradiction. A substantial meridional flow and re- 
distribution of magnetic flux is predicted at 4 AU, contrary to 
observation. Apparently the main reason for this is that the pres- 
sure in the solar wind is anisotropic and, together with a consid- 
eration of how the effective pressure may vary with heliocentric 
distance, we find that Ulysses and other data for the pressure 
anisotropy are entirely consistent with only a small meridional 
redistribution of magnetic flux at 4 AU. This is, therefore, an 
example of a breakdown in the simplest form of the isotropic 
fluid equations in application to the large scale dynamics of the 
solar wind. 

2. Analytic estimate of meridional transport due to a 
meridional thermal pressure gradient 

We consider the solar wind to be described by the ideal MHD 
equations for a single fluid and we apply these to analyzing the 
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extent a merdional thermal pressure gradient can redistribute 
magnetic (and mass) flux in the interplanetary medium between 
10 Rq and 1-4 AU, a region in which 3 > 1. Assumptions 
can be made for this specific problem which make the analytic 
treatment tractable and simple. Including those stated above, 
these are: 

(1 ) The flow is nearly (but not exactly) radial. 

(2) The radial flow speed is independent of radius. 

(3) Magnetohydrodynamic effects are negligible. 

(4) Gravity is negligible. 

(5) The flow is axisymmetric. 

Assumption (1) is valid beyond 5 Rr. ly except during large 
transients (which are not under consideration here). (2) is suf- 
ficiently well satisfied to permit an estimate of whether merid- 
ional mass and magnetic flux transport is a large effect at 1-4 
AU. It is not necessarily sufficiently accurate to permit a quanti- 
tative calculation of meridional transport. (3) has been evaluated 
elsewhere to show that MHD effects are smaller than, and hence 
can be treated separately from, thermal pressure gradient effects 
(Suess and Nerney, 1975a,b). (4) follows simply by knowing the 
solar wind flow is supersonic and from comparison of the ki- 
netic energy with the gravitational potential energy. With regard 
to (5), the Ulysses fast latitude scan in 1994-1995 was just prior 
to solar minimum, when the Sun’s magnetic field might rea- 
sonably be approximated as an axisymmetric dipole. At other 
times, the results we find here would only be applicable, at best, 
to properties averaged over a solar rotation period. This may be 
a weak assumption when there exist large corotating interaction 
regions, as shown by Pizzo and Goldstein (1987). 

Under the above assumptions, the solution to meridional 
transport can be estimated using a “quasi-radial flow” approxi- 
mation (Suess and Nerney, 1975a). Consider B, v, «, and T to 
be written as 


B = Bi -h B 2 + ... (1 a) 

v = \\ + v 2 + ... (\b) 

n — n\ + n 2 4 - ... (lc) 

T — T\ + Tj 4 - ... (!</) 


where vi = e r v\ r and it is assumed that | B? |<C| B] |, 
I V 2 |<| vi |, etc. Terms of 0[1], those with subscript ‘ 1\ 
are collected together to produce Eqs. (2) and the 0[1] solu- 
tion. Terms of 0[B 2 , v 2 .rt 2 , 73] are collected to find the next 
“higher order” equations. The derivation shown here is fully and 
completely described below, but does not include some “back 
pressure” terms which would limit flux redistribution in lati- 
tude if it becomes an 0[1] effect. Such might be the case in the 
outer heliosphere and will be the subject of a more comprehen- 
sive analysis elsewhere. A general discussion of perturbation 
expansions of differential equations is given, for example, by 
Cole (1968). 

The “lowest order” equations describing the solar wind in 
the present approximation are then c ]r = V r (0) and 

-t-j -(r 2 n\(r.0)V r {0)) = 0 (2a) 

dr 


1 d , 

— — (#-*, r (r)) = 0 (2b) 

r L dr 

or 

ni(r.O) = n E (0) (y) (2c) 

B\ r (r) — B or (y) (2 d) 

where B or is the value at the top of the corona, r (} =10 Rq 9 
beyond which distance 3 is greater than unity, and the density 
at 1 AU, n E (0), is from Ulysses observations scaled to 1 AU= 
r E . Because of the results from Suess and Smith (1996), it can 
be taken that B or / a function of 0. The parameters defining 
n E (0) and the temperature will be determined empirically, using 
Ulysses data. There is no momentum equation to this order 
because V r (0) is independent of radius. There is a meridional 
pressure gradient, 8p\j80 , but its dynamic effect in producing 
a meridional flow is shown below to be scaled by A/~ 2 , where 
M is the Mach number. Therefore, the meridional momentum 
equation is of the next higher order, and is used to determine the 
the absolute magnitude of the higher order effects. The resulting 
meridional flow in turn modifies the mass and magnetic flux 
distribution. 

Four equations will be involved in computing the merid- 
ional transport - the meridional component of the momentum 
equation, the continuity equation, the solenoidal condition, and 
the induction equation. Eqs. (2) described the lowest order so- 
lution. In the next order, the linearized meridional momentum 
equation is given by: 

m p ni(r,e)^-^~(rv 2 g(r,e)) = (3) 

r or r 80 

where the formula has been written in such a way as to em- 
phasize that the lowest order flow field (subscript ‘ V) provides 
the forcing function to the next higher order (subscript ‘2’) re- 
sponse of the solar wind. This ordering only works if the so- 
lution CT 2 , B 2 , v 2 , n 2 ) remains small relative to the lowest order 
solution. The rest of the equations in this order are 

vioir. 0)B] r (r) - V r (0)B 2e (r , 0) = 0 (4a) 

— «-( r 2 B 2r (r , 9)) + — ]— ( B 2 g(r , 6) sin 6) = 0 (4b) 

r 1 or r sin 6 86 

4 S-(r 2 n 2 (r. <9)V r (6>»+ 
r" or 

1 8 

— o(r,0)s\n0) = 0 (4c) 

r sin 0 80 

Eq. (3) can be seen to be a valid linearization in the su- 
personic solar wind if it is written in dimensionless form by 
scaling the velocity, density, and the temperature with their val- 
ues (U] r , «£> Te) at 1 AU for an arbitrary value of 0 — 0 E . Then, 
with primes denoting dimensionless variables 

d(r'v f 29 (r\0)) _ 1 1 8p\(r\0) 

dr' ~ M 2 E p'\> r (0) 00 


( 5 ) 
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In words, the 0[1] merdional gradient in p] leads to an 
0[1 jM\] meridional flow, v' 10 (r f ,6), where M\ is the square 
of the Mach number at 1 AU, at the polar angle 0 = 0 £ . Since 

>> 1 in the supersonic solar wind, [v' 20 (r\ 0)] << [V'(0)]. 
There are, of course, other terms in the full meridional momen- 
tum equation. These include a term of 0[ #)], which is 

of 0[1 jM\\ and hence ignorable. They also include a term of 
order which would have to be as large as Ml to be 

important. This is the back pressure referred to earlier which 
limits the redistribution of flux w'hen it becomes large. Other 
terms are ignored under the assumptions listed at the beginning 
of Sect. 2. With regard to the rest of the equations which deter- 
mine V 2 o ( r < 9), (4a) is the azimuthal component of the induction 
equation, and (4b) and (4c) are the solcnoidal condition and the 
continuity equation, respectively, to this order. 

The solution procedure is the following. Empirical functions 
for V r (9), p\(r, 0 ), n\(r, 0), and T](r, 0) are used in (3) to com- 
pute V 2 e(r< 0). This is used in (4a) to find B2o(r<9). Eqs. (4b) 
and (4c) are then solved for «2(r,0) and #2r(A 8). The latter is 
the desired dependent variable since it describes the redistribu- 
tion of magnetic flux by the meridional pressure gradient. The 
overall result is given by 

B r (r.0) = B lr (r) + B 2 r(r,9) (6) 

It remains to specify the empirically determined analytic 
forms of p\{r, 0 ), n\(r , 6 ), V r (0), and 7j(r, 0). 

3. Application of Ulysses fast latitude scan observations to 
the evaluation of meridional transport 

We first turn to the empirical determination of n\ and 7j, and 
we also combine the plasma data with Ulysses magnetometer 
data to show that 3 > 1, justifying assumption (3) above. The 
Ulysses solar wind plasma instrument measures the ion and elec- 
tron plasma distributions, from which the temperature, density, 
and pressure can be determined (Phillips et al. 1995). Ion and 
electron temperature radial gradients are not well determined, 
but other empirical determinations give the non-adiabatic de- 
pendence of these temperatures on radius as approximately 
r -° 7 (Phillips et al., 1995). Here, we will adopt this depen- 
dence but will later discuss the effect of other choices such as 
the r~° 5 dependence found by Richardson et al. (1995) in the 
IMP 8 and Voyager 2 data. Any empirical radial dependence is 
easily utilized in the analytic evaluation. 

The total pressure of the plasma is then given by p — 
k(n p T p +n a T a ^n l T e ),whcrc(n p ,n a ,n e )2ind(T p , T a , Tjarethe 
number densities and temperatures of protons, alpha particles, 
and electrons, respectively. n e is, occasionally, not accurately 
measured and therefore is replaced by (n p + 2 n^) — n e . The 
proton and alpha particle temperatures are the one dimensional 
temperatures in the energy direction over full field of view. This 
is the most accurate temperature available from the Ulysses 
SWOOPS data. 



Fig. 3. 25.5 day running averages of the thermal pressure at Ulysses, 
scaled to 1 AU as described in the text. Also shown are the heliographic 
latitude of Ulysses and the curve fit to the pressure data that is described 
in the text. 

Using the above radial variation, the observed total ther- 
mal pressure can be scaled to 1 AU. This quantity is subject to 
large fluctuations due to transient disturbances. Therefore, for 
the purposes of studying global scale meridional redistribution 
of mass and magnetic flux, running averages over 25.5 days (one 
solar rotation) are used. Averages are computed as the average 
of products. The result is shown in Fig. 3, along with the heli- 
ographic latitude of Ulysses, for the interval of the Ulysses fast 
latitude scan between 1 August 1994 and 1 September 1995. 

The observations shown in Fig. 3 were made over a brief 
one year interval so that global solar conditions are relatively 
constant. Plasma data from this interval have been reported in 
more detail by Phillips et al. (1995). The normalized pressure 
can be seen to be of Of 1 .5 - 2] x 1 0~ 10 dynes/cm 2 over most 
of the interval, being approximately constant from 80.22° S to 
50° S and again from 50° N to 80.22° N. The pressure probably 
continues at about the same value all the way to the North and 
South poles because at the higher latitudes Ulysses was fully 
immersed in the polar high speed streams. However, between 
50° S and 50° N the pressure increases by about 50%, maximiz- 
ing at the equator. Therefore, there is a pressure maximum at 
the equator with a relatively narrow half-width that is approx- 
imately symmetric across the equator. A curve has been fit in 
Fig. 3 (dashed curve) to the pressure variation. Using a scaling 
in radius of density oc r~ 2 and temperature oc r~ 01 (Phillips 
et al., 1995) to define the radial variation of the pressure, the 
equation for the curve fit is given by 

P\(r.e)=p E (\ +e<r acc,s:s ) (y)’ 7 (7) 

The parameters for the fit shown in Fig. 3 are = 1 .44 x 1 0“ 10 
dyne-cm 2 , a = 3.5. and e — 0.9. This will serve for examining 
whether meridional pressure gradients might cause significant 
mass and magnetic flux redistribution. 

Before applying (7), we first consider the plasma 3, found by 
combining the pressure shown in Fig. 3 with the total magnetic 
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Fig. 4. 25.5 day running averages of the plasma 0 at Ulysses during the 
fast latitude scan, scaled to 1 AU. The radial magnetic field is scaled 
as r “ 2 and the transverse field as r~ 1 . The temperature and density are 
scaled as described in the text. 
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Fig, 5. Density and radial flow speed versus time for Ulysses during the 
fast latitude scan. The density is scaled to 1 AU as r -2 . 25.5 day running 
averages of the data has been used in both the density and velocity plots. 
Curve fits to the data are superimposed, with the parameters given in 
the text. 


field strength at Ulysses, scaled to 1 AU, with a running average 
over 25.5 days. These results are shown in Fig. 4, for the interval 
covered in Fig. 3. Over the entire interval ,3 > 3 and it can 
therefore reasonably be asserted that MHD effects are relatively 
small, as required by assumption (3). 

Now, carrying out the integration of Eqs. (3) and (4) de- 
scribed in Sect. 2 results in the following solution: 


At this point, we can now answer the posed question; how 
much does the observed meridional thermal pressure gradient 
redistribute magnetic flux between 10 Rq and 1-4 AU? 

To do this, it is easier to write (8c) in terms of the ratio 
between B 2r and B \ r , to give 


ggrOvg) 

B\ r (r) 


2nep E (r E /r„) 01 

O.3« £ (0)m„V r (0) 2 


V2o(r.9) = — - sin 2 6e QC0S e 

03m p n E (9)V r (9 ) 


* ( r 


03m p n E (6)V*(0) 


xf^-J (r°- 3 -r™) 


B 2r (r,9) = 


latpE^r^rf^ 

03n E (9)m p Vj(9)r 2 




x 2 cos 2 9 — sin 2 9 — 2a sin 2 9 cos 2 9 


1 2 ‘'McoJ,--'* 


n E d9 V r d9 


r \ 07 

~ 1 


x 2 cos 2 9 — sin 2 9 — 2a sin 2 6*cos 2 9 

4iS + ^) co 4'“‘“ : ’ <9) 

This result depends on the individual velocity and density merid- 
ional gradients, in addition to the overall meridional pressure 
gradient. The density and velocity variations are shown in Fig. 
5 in the same manner as for the pressure in Fig. 3. 

That is, 25.5 day running averages are made. The density is 
scaled to 1 AU as r~ 2 and the velocity is assumed independent 
of radius. Then the two sets of data are fit by the same functional 
form used for the pressure: 


V r (0) = v 0 ( 1 4- e v e 


n E (r , 9) = n oE ( 1 + t N e 


)S 2 

* v cos : 0 ^ 


The curve fils in Fig. 5 are for v a = 771 km/s, e v — -0.45, 
a v — 5.5 and n oE = 2.85 cm -3 , e n = 1.73, a„ = 5.5. These 
are substituted into (9), which can then be abbreviated by 


B 2r (r.9) 


= Cfi(r)f 2 (9) 
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Heliographic Latitude 


Fig. 6. The relative variation of the modeled magnetic field strength 
versus heliographic latitude at 0.25, 0.5, 1 .0, and 4.0 AU and at r — ► oc 
(from top to bottom). This is from Eq. (11). 


where C = 0.92, to give 



Fig. 7. One hour averages of the radial flow speed over the interval of 
the Ulysses fast latitude scan. The slot around the heliographic equator, 
which lies over the solar streamer belt, is seen here to have a series of 
high speed streams. 


B 2r {r.d) 

B\ r (r) 


0.28 f 2 (0) 


(12 a) 


B 2r (r.9) 
B\ r (r ) 


0.39/2(0) 


(12 b) 


The result suggested by Eq. (12) is significant. It is that flux 
redistribution due to the pressure, density, and velocity merid- 
ional gradients illustrated in Figs. 3 and 5 is not negligible. Fig. 
6 shows a plot of Eq. (9) (or, equivalently, Eq. (1 1)), for the 
given parameters, for the relative variation of the absolute field 
strength at 0.25, 0.5, 1 .0. and 4.0 AU and infinitely far from the 
Sun. A broad, deep deficit in radial magnetic flux is predicted 
around the equator, with a magnitude of approximately 25% far 
from the Sun. Not only is B 2r /B\ r not negligible, but also v 2 $ 
given by (8a), using the above parameters, is not negligible. The 
maximum value of v 2 o, which occurs about 35" away from the 
equator, is 17 km/s at 1 AU and 8 km/s at 4 AU. A meridional 
flow this large would probably be observable by the solar wind 
plasma instrument on Ulysses. 

We now must ask why the predicted result based on simple 
hydrodynamics differs from the observations. With regard to 
the linearization, M\ » I and v 2 e << v\ r throughout the 
solution domain, so the calculation is a mathematically valid 
estimate of the effect of the meridional pressure gradient. 

There are important reasons why the simple hydrodynamic 
result may not be applicable to the solar wind. One reason has to 
do with assumption (5), which was that the flow' is axisymmetric, 
or that stream interaction effects are negligible. Fig. 7 is a plot 
of the 1-hour solar wind speed measurements made at Ulysses 
throughout the fast latitude scan. A wide range of flow speeds is 
seen to occur in the region around the heliographic equator. This 
is precisely the same location as the pressure bump shown in 
Fig. 3. There are therefore corotating and/or transient high speed 
streams throughout the latitude range of the streamer belt. The 
equatorial measurements were made at around 1 .4 AU, where 


the streams are expected to have undergone some dynamic in- 
teraction and steepening. This is probably the source of some 
of the temperature (and hence pressure) increase in this range 
of heliographic latitudes. It was explicitly assumed above that 
the pressure bump originated at the Sun and exists continuously 
out into the solar wind. If, however, the pressure bump is partly 
due to stream interaction then the bump will not develop until a 
distance that is a large fraction of an astronomical unit. This pro- 
cess and the associated meridional transport has already been 
investigated in some detail by Pizzo and Goldstein (1987) in 
a numerical model of stream interaction. They found what has 
been shown again here analytically - that a narrow, steep-sided 
bump in the structure of the equatorial solar wind is able to cause 
a strong meridional transport of magnetic flux. Their calcula- 
tions began at 1 AU for the specific reason that stream inter- 
action is relatively negligible inside that distance. The Ulysses 
measurements during the fast latitude scan were made between 
1 .34 AU at perihelion, 2.29 AU at -80.22 degrees heliolatitude, 
and 2.02 AU at +80.22 degrees. It can be expected that stream 
interaction is therefore somewhat more important than at 1 AU. 
In terms of the present model calculation, it means a revision in 
the assumptions and the analysis would have to be made to take 
into account that the pressure bump had a different (smaller) 
relative magnitude closer to the Sun. 

To simulate the effect of an equatorial pressure bump which 
increases with heliocentric distance, the radial dependence of 
the temperature which enters into/i(r) in Eq. (11) can be 
changed. This can be used to examine the effect of heating due 
to stream interaction. In the present case, T — 2 x 10 6 K at 10 
Rq . The problem with this is that this temperature is probably 
higher than exists at that height in the corona because the extrap- 
olation was based on an in situ temperature that included some 
heating due to stream interaction. Therefore, the effect of the 
meridional pressure gradient near the Sun is overemphasized. 
The pressure given by (7) therefore overemphasizes the flux 





310 


S.T. Suess et al.: Latitudinal dependence of the radial IMF component - interplanetary imprint 




Fig. 8. The assumed proton temperature anisotropy at 1-4 AU as a 
function of latitude during the Ulysses fast latitude scan. The anisotropy 
reaches 1.5 at the equator and falls to 1.1 at the poles. The same variation 
is assumed to apply to alpha particles. The electron anisotropy at the 
time, also shown, is taken to be 1 .3 at all latitudes. These assumed 
values are used to compute the effective perpendicular pressure from 
the observed pressure shown in Fig. 3. 

redistribution effect near the Sun. Simply changing the expo- 
nent for the temperature dependence can partially simulate this 
effect. In the extreme case, the temperature could be assumed 
independent of radius. In this case it turns out that in Eq. (11), 
/i (r) oc ln(r)! This means that while there is a relatively smaller 
flux redistribution inside 4 AU, there may be a large effect in 
the outer heliosphere. An alternative assumption might be to 
take half the pressure bump to be due to stream interaction. The 
result in this case is that the effect in Eq. (1 1) is simply halved 
as well. 

Therefore, stream interaction can reduce the flux redistribu- 
tion, but does not seem to be able to reduce it to the observed 
levels for strong corotating interaction region heating. Further- 
more, it is not obvious that corotating interaction regions are 
producing much heating in CIRs during the the Ulysses fast lat- 
itude scan. The streams shown in Fig. 7 are small in comparison 
to the model streams invoked by Pizzo and Goldstein (1987) to 
illustrate CIR effects. Conversely, streamers are known to pro- 
duce high temperatures and densities in the HPS, as shown by 
empirical results and streamer models (J. Phillips, priv. comm., 
Suess et al., 1996). 

What apparently has a far more important effect on the pres- 
sure gradient driving meridional flux redistribution is the pres- 
sure anisotropy. As noted above, the proton and alpha particle 
temperatures used in computing the pressure shown in Fig. 3 
are the one dimensional temperatures in the energy direction 
(approximately the radial direction) over the full field of view. 
However, it is the perpendicular temperature in the meridional 
direction that gives the effective meridional pressure gradient. 
The electron, proton, and alpha particle anisotropies are mea- 
sured by SWOOPS on Ulysses and, although the data have not 
yet been fully reduced, preliminary results indicate that the elec- 


Fig.9. The 25.5 day running average perpendicular proton temperature 
during the Ulysses fast latitude scan using the assumed anisotropy 
shown in Fig. 8. The thin line is a curve fit to the data using Eq. (2), with 
the parameters q = 3.33, e = .58, and pE = 1.25 x 10~ ,() dynes/cm 2 . 
The isotropic pressure from Fig. 3 is also shown for reference. 

tron anisotropy is approximately 1 .3 everywhere, while the pro- 
ton anisotropy is approximately 1.1 in the high speed flow and 
1.5 in the low speed flow (Goldstein et al., 1995; Goldstein, 
Neugebauer, and Smith, 1995; Goldstein et al., 1996). 

An extensive previous study of proton anisotropies was re- 
ported by Marsch et al. (1982) using Helios 1 and 2 data. In 
Fig. 15 from that paper, it can be seen that the Ulysses results 
are consistent with the Helios results as extrapolated to 1 -4 AU. 
Inside 1 AU, they found that the anisotropy in the high speed 
solar wind was often less than unity, while it was as large as 2.0 
in the low speed flow. These trends will be used below to show 
that using the values in Fig. 8 , independent of radius, gives an 
upper bound to the amount of flux redistribution. 

The analytic estimates in Eqs. (3) and (4) were made using 
an isotropic temperature and pressure. Because the temperature, 
and hence pressure anisotropy is larger than unity, the perpen- 
dicular pressure will be reduced relative to that shown in Fig. 3. 
Beyond this simple effect, the difference in pressure anisotropy 
between the high speed and the low speed solar wind means that 
the pressure gradient, or difference, between these two regimes 
will be further reduced. Hence, the effective meridional pres- 
sure variation is considerably smaller than shown by Fig. 3. 
Applying the anistoropies in Fig. 8 to the electron, proton, and 
alpha particle pressures used in Fig. 3 gives the result shown 
in Fig. 9. 

This figure shows that the effective pressure for computing 
the meridional pressure gradient is far smaller than would be 
expected by considering only the isotropic pressure. Using again 
Eq. (7) to approximate the pressure, we now find that a = 3.33, 
e = 0.58, and p E = 1 .25 x 10“ 10 . The curve derived from these 
parameters is shown as the thin line in Fig. 9. Substituting these 
revised numbers into (11) gives C = 0.49 and 

~ p r 7 . ~ \\au^ OA5f 2 (0) (13 a) 
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Hellographic Latitude 


Fig. 10. The relative variatoin of the modeled magnetic field strength 
versus heliographic latitude at 1.0 AU for the perpendicular pressure 
shown in Fig. 9 (solid line). The relative variation using the isotropic 
pressure, also at 1 AU and taken from Fig. 6, is shown by the dahed 
line for comparison, t = .58, and pf. = 1 -25 x 10 10 dynes/crrr. The 
isotropic pressure from Fig. 3 is also shown for reference. 


* 2 r(r,fl) 

B] r (r) 


0 . 21 f 2 (0) 


(13 b) 


for the flux redistribution using the perpendicular pressure 
shown in Fig. 9. This is a much smaller effect than found in 
Eq. (12). A comparison of the calculated flux redistribution at 
1 AU using the pressure shown in Fig. 3 with that found using 
the perpendicular pressure in Fig. 9 is shown in Fig. 10. The 
comparison is made at 1 AU and it is seen that the assumed 
anistropies shown in Fig. 8 result in about a 50% reduction in 
the computed flux redistribution. 

Similarly, the maximum meridional flow speed is also re- 
duced by about 50%, from 17 km/s to 9 km/s at 1 AU, at a 
latitude of 35". 

The pressure anisotropy used here is independent of radius. 
If the radial variation observed by Marsch et al. (1982) had 
been invoked, the flux redistribution would have been reduced 
by approximately another 50% . This is because the difference 
in the perpendicular pressure between the high speed and low 
speed solar wind is smaller closer to the Sun than it is at 1-4 
AU. This implies that the pressure anisotropy alone is enough 
to account for the discrepancy between the original calculation 
shown in Fig. 6 and the observations of IMP, ICE, and Ulysses. 


4. Discussion 

The conclusion from the analysis in Sect. 2 and Sect. 3 is that the 
meridional pressure distribution observed by Ulysses at 1.4-4 
AU will not modify the distribution of magnetic flux already 
present in the solar wind at the top of the corona - 10 Rq . 
This conclusion would not have been possible if the pressure 
anisotropy were ignored. 

With regard to the consistency of the analysis, reconsider the 
assumptions listed in Sect. 2. All of them seemed well-justified 
except (2) - that the radial flow speed be independent of radius. 


This approximation depends on how much solar w ind acceler- 
ation occurs inside 10 R . Recent studies (e.g. Habbal et al 
1995) have suggested most of the solar w ind wind acceleration 
is complete by this distance so the approximation is probably 
good. Classical solar wind models without added acceleration 
source terms have a much flatter speed profile in radius so that 
approximation (2) would be much less valid. However, in the 
end, the conclusion here that the meridional pressure gradient 
does not lead to a large magnetic flux redistribution means that 
the details of the flow speed profile in radius are not critical to 
the conclusion. 

A critical assumption could have been ignoring stream inter 
action. This is not because of the driving of meridional flow due 
to pressure ridges associated with passing high speed streams 
That process has been examined by Pizzo and Goldstein ( 1 987 ) 
and found to be negligible. Instead, the heating by stream in 
teraction can be an important source for the pressure bump 
This means that near the Sun the pressure bump might have 
been smaller, and therefore the effective meridional transport 
reduced. However, the Ulysses fast latitude scan data was taken 
near solar minimum when the amount of heating by CIRs does 
not appear to be very large. Furthermore, coronal models imply 
that there will be a significant pressure bump in the signature of 
the streamers themselves (Suess et al. 1996). The dependence 
of/i(r) in (1 1) on radius can be used to investigate this, other 
limiting cases, and the consistency of the analysis. Taking the 
extreme case of a temperature being independent of radius gives 
f](r) oc ln(r) and there is a secular increase of Bi, with radius. 
This case turns out to be analogous to redistribution due to inter 
nal MHD forces, where it was shown earlier (Suess and Nerney. 
1973, 1975b) that the MHD effect is proportional to ln(r). In the 
MHD case meridional flow is driven by the gradient in B 2 a with 
latitude. Asymptotically, B 0 cx (1 /r) and the equivalent mag- 
netic pressure is proportional to (1 /r) 2 . This is the same as the 
dependence of the thermal pressure on radius for an isothermal 
wind. Therefore, the two results are consistent with each other. 
The magnitude of the MHD effect is nevertheless very small in 
the inner solar system for the simple reason that S3 > 1 through 
out the region under discussion. This places an upper bound, 
also, on what might be expected from a weaker radial depen 
dence of temperature on radius, such as found by Richardson et 
al. (1995). 

The pressure anisotropy is found to be the most importani 
physical parameter for understanding the Ulysses observation^ 
of the distribution of B r in latitude. The perpendicular pressure 
is smaller than the isotropic or parallel pressure, as suggested 
in Fig. 9. It is the perpendicular pressure that drives meridional 
transport and the anisotropy reduces the effective transport to 
about 50% of that predicted using an isotropic pressure. Con 
sidering that the pressure anisotropy shown in Fig. 9 is also a 
f unction of heliocentric distance and using the empirical results 
for this radial dependence given by Marsch et al. (1982) shows 
that meridional transport near solar minimum can be expected to 
be < 10% at 1 AU and < 20% at larger distances. It was shown 
earlier (Suess and Smith, 1996) that the field strength at 10 7? 
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is constant except in the HPS so D r should be independent of 
latitude on global scales between 10 R ( - } and a few AU. 

With regard to the quasi-radial How approximation itself, 
recall that Eqs. (2), (3), and (4) have been linearized around 
radial flow and the meridional gradients of p 2 and B 2r have 
been ignored in computing v 2 o- If these gradients become too 
large, they will resist continuing meridional flux redistribution. 
It can be, in other words, ultimately a self-limiting process if it 
continues rapidly in comparison to a solar wind expansion time. 
The amount of flux redistribution computed here during the 
Ulysses fast latitude scan is insufficient to lead to a significant 
back pressure at any location inside 10-20 AU so that this is not 
an important effect for the present study. 

It is possible to reconstruct the analysis taking the back- 
pressure force into account, but this would have added unnec- 
essary complexity to addressing the fundamental question of 
whether the problem might produce an interesting result in the 
first place. It is still possible that there will be significant merid- 
ional mass transport in the outer heliosphere due to the separate 
processes of heating due to corotating interaction regions and 
the meridional dependence of the IMF. However, that is not the 
topic of the present investigation. 
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